multipleOLS <- function(formulas,y=data,modelnames,...){
	require(lmtest)
	foo <- vector("list",length(formulas))
	for(i in 1:length(formulas)){
		foo[[i]] <- lm(substitute(i,list(i=formulas[i])), data=y)
		foo[[i]] <- coeftest(foo[[i]],vcov=vcovHC(foo[[i]],"HC1"))[2:nrow(coeftest(foo[[i]])),1:2]
		foo[[i]] <- cbind(foo[[i]],
 #           foo[[i]][,1]/foo[[i]][,2],
						foo[[i]][,1]-1.96*foo[[i]][,2],
						foo[[i]][,1]+1.96*foo[[i]][,2])
		foo[[i]] <- data.frame(foo[[i]],modelnames[i],rownames(foo[[i]]))
		colnames(foo[[i]]) <- c("Estimate","SE","Y0","Y1","Model","Var")
			}
	return(foo)
}
